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r— I The Tevatron heavy flavor physics program is in fuU swing. The rapid increase in the size of 

data samples is allowing significant improvements of previous results, and opens the doors to 
new possibilities. A further doubling of the current integrated luminosity is expected in the next 
couple of years. This report summarizes the main current results and future prospects. 
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1. Properties of Heavy-Flavored hadrons 

The rich production of heavy hadrons of all kinds, together with limited event pile-up and 
high-quality tracking detectors puts the Tevatron in a very favorable position for many experimen- 
tal studies of heavy-flavored hadrons. Some topics are even exclusive to Tevatron. One example 
in this category is the properties of bottom-flavored baryons. Just few years ago, the was the 
only bottom baryon for which some significant experimental information was available. In the past 
couple of years, knowledge of Z7-baryons has expanded enormously thanks to Tevatron data. In year 
2007, CDF observed the strongly-decaying {uub) and {ddb) baryons\ together with their 
orbitally excited (/ = 3/2) partners Z^^, Z^~, and measured their masses with a precision at the 
MeV level[3]. That same year, both D0 and CDF observed the charged baryon {usb), making 
the not anymore the only known example of weakly-decaying bottom baryon [4, 5]. Finally, 
last year D0 announced the observation of the doubly-strange Q.^ baryon {ssb), the heaviest fun- 
damental state of a single-bottomed baryon, recently followed by a similar result from CDF. The 
masses of heavy-flavored baryons are more than just a "zoological" curiosity: they allow to test 
our understanding of the dynamics involved in the structure of hadrons. The availability of a whole 
spectrum of measurements from the Tevatron is particularly important in this respect, as most mod- 
els are much better at predicting differences of masses or other overall relationships than individual 
masses. Here we discuss the most recent results, that are about the weakly-decaying baryons 
and Q.^. Both states are easily triggered and reconstructed at both CDF and DO through their de- 
cays into J/\if IJi^lJi~ : J/\ifQr , Qr APk~ and J/W'^~ ' 2" ^ A^tt". The 
presence of multiple long-lived particles makes for a three-vertex, 5-track configuration. In the 
case of DO, a special-purpose tracking code has been used to improve efficiency for reconstruction 
of tracks with significant impact parameters to the beamline. 

The initial measurements of the baryon mass are shown in Fig. 1. They are in good agree- 
ment with each other and with predictions, although the CDF measurement is significantly more 
precise, matching the resolution of the most precise predictions available. CDF has recently per- 
formed an update with a 4.2fb~^ sample, performing a global analysis aimed at both the 3^ and Q.^ 
baryons, and including several other B meson channels of similar decay topology as controls [6]. 
The result: M(H^) = 5790.9 ± 2.6 =b 0.8MeV/c^ is compatible with the previous measurement and 
supersedes it. For the CDF finds M(t2^) = 6054.4 ± 6.8 ±0.9MeV/c^. This level of precision 
allows to actually discriminate between different models: for instance, the new masses measured 
by CDF are in close agreement with recent lattice [8] or color-hyperfine-splitting calculations [9], 
with a resolution that appears to allow even a discrimination between different forms of the po- 
tential. However, there is a 6(7 disagreement between the CDF measurement of M(f2^) and the 
value previously measured by DO: M(f2^) = 6165 ± 10 ± 13MeV/c^. In addition, the production 
rate is also not very well matched, although both experiments claim signal significances above 5(7, 
(see [6, 7] for more details). This puzzle must be solved before definitive conclusions can be drawn, 
and this will likely not happen without the help of further experimental data. 

Both CDF and DO find the decay time distribution for the to be compatible with a weakly- 
decaying Z7-hadron. CDF also quotes explicit lifetime measurements. The lifetime is measured 
exclusively for the first time, yielding: t(S^) = 1.56^q 25 ps; the result is consistent with 

^N.B. the two oppositely-charged states are not a particle-antiparticle pair. 
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Figure 1: masses (left) and Q.^ signals (right) 



the average lifetime obtained at LEP for a mix of S^, S^, and baryons [1]. For the Q.^ this 
is the first lifetime measurement: T(fl^) = 1.13^q 4q ±0.02 ps, again compatible with a weakly- 
decaying b hadron, and slightly favoring a shorter lifetime than the B^. In addition to these unique 
measurements on bottom baryons, many other measurements of bottom hadron lifetimes have been 
performed at the Tevatron. The lifetime of the unique doubly heavy-flavored meson Be has been 
measured with 10% resolution [10], and the current world averages of and 5^ lifetimes are 
dominated by the Tevatron measurements performed with lfb~^ [34, 2]. All such measurements 
are expected to improve with the large samples now available and not yet analyzed (^ 5fb~^). Even 
the B^ and B^ lifetimes, for which large samples are already available from e^e~ experiments, will 
be further improved and become dominated by the Tevatron measurements. 

Progress has also been made in the field of exotic heavy hadrons. CDF has the largest exist- 
ing sample of X(3872), the first and most studied exotic state, allowing precision measurements 
aimed at understanding the nature of this still mysterious particle. In the past these have allowed to 
strongly constrain the possible J^^ assignments [12]. The latest analysis of a 2.4fb~^ sample, cor- 
responding to c:^ 6,000 signal events, yields the most precise mass measurement: M(X(3872)) = 
3871.61 ±0.16it 0.19MeV/c^, and allows a detailed analysis of its lineshape, excluding the possi- 
bility of a two-component system with a separation greater than 3.6 MeV/c^ in the assumption of 
equal production [13]. The measured X(3872) mass is lower than M{D^) +M(D*^), still allowing 
the possibility of interpretation as a loosely bound "molecular" system. However, the masses are 
compatible within uncertainty, the difference being just 0.19 =b0.43MeV/c^, with the uncertainty 
being dominated (somewhat paradoxically) by the PDG uncertainty on the masses of the common 
mesons. The above observations favor molecular models over tetraquark models; however, the 
observed large production cross section of the X(3872) at the Tevatron seems very difficult to rec- 
oncile with a molecular model [14], so the issue is far from being settled, and the true nature of the 
X(3872) remains a fascinating puzzle. 

The study of further states may help understanding what is going on. The latest appearance 
on the scene is the Y(4140), a 3.8(7 excess observed by CDF in the mass spectrum with 

mass M(F(4140) = 4143.0 ± 2.9 ± 1.2MeV/c^ [15]. This has been reconstructed in a sample of 
15B^ J/\l/(j)K^ decays (also the largest existing sample of this mode) obtained from a 2.7fb~^ 
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Figure 2: Current results (left) and constraints imposed by some measurements on the MSSM model (heavy 
squarks) [23, 24] (right) 

sample. The nature of this signal is still unclear and will benefit from the addition of further data. 
2. Rare Decays 

The large B and D production rate and the strong background rejection capabilities of the 
two detectors make the Tevatron an ideal place to look for rare heavy flavor decays. Several rare 
modes have been studied, but the channels attracting the greatest interest are 5^ and 
. These Flavor Changing Neutral Current processes have small branching fractions in 
the Standard Model, but may receive contributions from a large variety of BSM processes. The 
latest and most precise predictions of branching fractions in the Standard Model are BR{B^ 

= (3.6 ±0.3) X 10-^ mdBR{B^ ^ i^^i^-) = (1.1 ±0.1) x 10-^^ [19]. It is worth noting 
that the ratio between the B^ and the B^ modes is not necessarily preserved in SUSY models, and 
the B^ might even be larger than the B^ (see for instance [20]), and it is therefore important to study 
both separately. 

Both CDF and D0 have now analyzed significantly larger samples since their last publica- 
tions [22], but the final results are still "blind". The new D0 analysis uses about 4fb~^ , divided 
into three different run ranges, that are analyzed independently, and then combined to produce a 
single limit with the CLs algorithm [38]. A Boosted Decision Tree technique is used to optimize 
the selection, by minimizing the expected upper limit under the assumption of no signal, resulting 
in an expected limit of 4.3(5.3) x 10"^ at 90(95)% C.L. 

The new CDF analysis is based on 3.7fb~^ and makes use of an Artificial Neural Network 
(ANN, or NN) to optimize the selection. Data are not divided into periods, but are instead divided 
in bins of NN output, to exploit differences in signal/background ratio to increase sensitivity. The 
new analysis has an increased geometrical acceptance with respect to the past, owing to the addition 
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of some further detector regions. The optimization procedure is otherwise the same as done by D0. 
The new expected limit for the 5^ is 4.3 x 10~^, almost a factor 2 better than previous 

results^. 

It is interesting to note that in this update CDF and D0 have an expectation for about one 
SM event each in their analyzed samples, although the background conditions prevent detecting 
it. The two experiments are now really closing in towards the Standard Model expectation, cutting 
into the final factor of ten, and will have significant impact in constraining the possibilities for 
the existence of new physics (see Fig. 2 for an example). Work is ongoing to include all data up 
to 5fb~^ and implement some further improvements to push the sensitivity of the analysis even 
further. A Tevatron average is likely to be very interesting at this stage. 

Several other results on rare modes are being produced in parallel with ji^ji' and 

Using data from the impact parameter trigger rather than the dimuon trigger, CDF has 
recently published results from the search of B^^B^ e^e~ ,e^ii~ that set the tightest existing 
limits on their branching fractions, and excluded the existence of leptoquark states with masses 
below ^50 TeV [17]. From the same trigger selection, tight limits on charm FCNC decays are 
also obtained [18]. Rare modes of the type B Xsi^i~ have already been observed at CDF 
with lfb~^ [25], and an updated analysis with the full statistics is ongoing^. These channels are 
particularly interesting in view of the recent hints for anomalies in the Aps distributions seen in 
e^e~ experiments [21, 37]. A confirmation of these anomalies might signal the presence of a 4th 
generation or other BSM phenomena. 

3. mixing parameters 

The B^ — 5^ system is an excellent lab for testing our understanding of SM physics and look 
for what may lie beyond that. Most conceivable new physics models have some impact on B^ oscil- 
lations, and the measurements of oscillation parameters work coherently together with the FCNC 
branching fraction measurements discussed in previous section to provide powerful constraints on 
BSM physics (see for instance [32, 19] and references therein). 

Experimentally observable quantities are: 

An, M M r + Ar r r r{B,{t) f) -m{t) ^ f) 

2 ' ' ^ r(B,(0^/)+r(B,(r)^/) 

They can be conveniently parametrized as follows: 

Ams = Amf^\As\ (3.1) 
Ar, = Arf ^ cos{<l)f^ + (j)^) (3.2) 

Ar, fAr,\'''cosi<t>!'' + ({>li ^^^^ 



Am. V Am. J A 



a%^^^t^n{<pf'' + <^^) (3.4) 

^CDF has since then unbUnded their result, obtaining: BR{B^ < 4.3 x 10"^ and BR{B^ < 

7.6 X 10~^ at 95% C.L. The expected Umit for B^ went down to 3.3 x 10~^, and the shght excess over 

prediction is quantified as a 0.73(7 effect [16]. 

^The analysis has been released at HCP 2009 and includes the first observation of the B^ (pji^ji' mode. 
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where |A^| and (^^ parameterize modulus and phase of possible non-SM contributions, and 
0™ is small (sin(0^^^) 2^ 0.004). In addition, we use 



(In the following, will omit the "// label for better readability). 

The current knowledge of the 5^ mixing parameters is dominated by the Tevatron measure- 
ments, and will stay that way until the LHCb will have collected and analyzed a sizeable sam- 
ple. The best known parameter is the oscillation frequency: Am^ = 17.77 ±0. 10 ±0.07ps~^ The 
experimental uncertainty is below the level of the theory uncertainty, and at the moment there 
is no obvious motivation for further improving it. The width difference is measured with lim- 
ited precision: AF^ = 0.062^q ps~^ [1]. This value is compatible with the latest predictions: 
AF^ = 0.088 zb0.017ps~^ [30], but it is also compatible with zero. However, this world-average 
was derived under the specific assumption of no CPV in the 5^ mixing, which is related to the 
mixing phase 0^. For this reason, and because of Eq. 3.2, a simultaneous determination of AF^ and 
0^ = —2^s is a better option. The phase 0^ is currently known only thanks Tevatron measurements 
of the time and angular distribution of ^ J/xf/ij), and still has a pretty large uncertainty. The 
asymmetry afs is hard to measure accurately. The only existing data is from the recent measure- 
ment by the D0 collaboration using semileptonic decays (A^^), yielding the result [31]: 

The DO capability of regularly inverting the B field direction contributes to keep the systematic 
uncertainty low. Although the best existing measurement, this is still compatible with zero, and far 
from testing the SM prediction of 2 x 10~^. 

The measurement of (AF^, 0s) parameters is based on an unbinned fit of the decay time and 
angular distribution of the J/ 1//0 decay. Both CDF and D0 have performed this measurement 
with 2.8fb~^; the reader is referred to other contributions in this conference for details [36, 35]. 
Here we report for the first time the Tevatron combination of (AF^, 0s) . This combination was 
performed by the Tevatron B averaging Working Group and includes CDF and D0 results with 
2.8fb~^ of data [26]. In order to produce this average, the CDF and D0 collaborations have worked 
together to adapt their initial measurements to make them compatible, making their procedures 
much more similar, to ensure the combined result is meaningful. As a consequence, the D0 result 
has changed slightly [35]. 

Both experiments now use a common method for the extraction of confidence contours in 
the (AF^, 0s) plane. Due to the critical impact that this measurement may have in revealing new 
physics, special care has been taken in avoiding approximations and staying on the conservative 
side if necessary. A complete frequentist construction with exact inclusion of systematic uncertain- 
ties is used. A profile likelihood function Lp^^/(AF^, 0s) is built for the two parameters of interest, 
by maximizing the likelihood over all other unknown parameters (physical parameters, like strong 
phases, and systematic parameters), and the ratio is taken with the value of the function at the 
mimum LRprof. The use of external constraints in the likelihood is minimized, and made equal 
between the two experiments. Due to limited statistics and strong non-gaussian effects, the usual 
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CDF Run II Prel. L = 2.8 fb ^ D0 Run II L = 2.8 fb" 




20log(/.) 20log(/.) 

Figure 3: Observed distributions of the profile HkeHhood ratio quantity for (AF^, 0s) measurement, com- 
pared to the nominal asymptotic chi-square distribution. 



assumption that the probability distribution of LRprof has the asymptotical chi-squared distribution 
is not a good approximation. The LRprof is therefore used just as the ordering function in a full 
Neyman construction of a multidimensional confidence region in the complete parameter space, 
which is then projected over the 2D space of the parameters of interests (AF^, 0s) . This procedure 
has been proved to produce confidence region with close-to-optimal perfomance in terms of sensi- 
tivity [27], and is critical to avoid an unnecessary weakening of the limits when the large parameter 
space (2^ 30 dimensions) is projected over the 2D space. The projection is performed by Monte 
Carlo generation of LRprof distributions for a number of points randomly sampled within the nui- 
sance parameter space, and choosing the distribution with the longest tail ("worst-case") as the 
basis for the construction of the confidence region. It is interesting to note that the deviation of the 
distributions obtained numerically from the nominal behavior is significant and cannot be ignored 
(Fig. 3). In order to avoid divergence of the size of the confidence regions, the allowed range of 
variation of the nuisance parameters was loosely bounded without compromising the coverage of 
the procedure, by exploiting the method of ref. [28]. Finally, the 2D confidence regions obtained 
in this way are combined with the usual HFAG method [2]. 

The final combined result is shown in Fig. 4. A 1-dimensional confidence interval for is 
also obtained: 



ft G [0.27,0.59] U [0.97, 1.30]@68%CL (3.7) 
ft G [0.10, 1.42] @95%CL (3.8) 

(3.9) 

The above confidence intervals imply that the SM value ft = 2.1 dz 0.7 x 10~^ is excluded at 
more than 95% CL. More precisely, a goodness-of-fit test of our data with the SM model , yields a 
significance of 0.034 (given the SM uncertainty on the value of Ay, the p- value for the "most con- 
servative" SM point is actually 2.0(7). Comparing this result with the average performed by HFAG 
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2008 reveals that the level of agreement with SM is essentially unchanged. This comes from the 
opposite effects of two changes: the increase of statistics on the CDF side, and a better accounting 
for tails in the D0 measurement. Further work is in progress to allow a unified simultaneous fit 
of the CDF and D0 data samples, which should in principle provide a better resolution because of 
a more effective use of the information in the parameters that get "profiled away" in the step of 
taking 2D projection of each experiment. 

In conclusion, the value of keeps showing a mild disagreement with the SM expectation. 
The fit prefers rather large values of j8^, ^0.5. The jury is still out on whether this is first sign of 
new physics (many BSM models predict large values of j8^, most notably 4th-generation models) or 
just a statistical fluctuation. The answer should become clear by the end of run II. It is interesting 
to note that in some models, a large phase implies a 5^ ^ /i+/i~rate just around the corner from 
current limits. 

3.1 5^ ^ 00 

Amongst the most recent results is the high-statistics study of the ^ (p(j) mode. This is a 
b sss pure penguin process and it is sensitive to possible new physics CP violating contributions, 
both in the decay and in the mixing. In addition to Acp, its BR and polarization are also good 
measurements to search for deviations from SM. 

This mode was first observed at CDF in 2005, and still unique to this experiment. The first 
observation was based on 0.2fb~^ and only had 8 signal events; the current analysis has exploited 
not only an increased sample of 2.9fb~^ but also an improved acceptance thanks to a better usage 
of multiple trigger selection, resulting in a sample of about 300 events (Fig. 5). This allows now a 
much more accurate determination of BR: 

Ks VV) ^ (i.78±0.14±0.20) X 10-2 ^.g^Q^ 



BR{BO,^J/Y(l>) 

BR{B^,^^^) = (2.40 ±0.21 ±0.27 ±0.82) x 10"^ (3.11) 
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Figure 5: B^^ ^ (jxj) signal at CDF 



where the errors are respectively due to statistics, systematics, and branching fraction uncer- 
tainty from the PDG. This compares well with predictions of a recent QCDF calculation: (2.18 ± 
0.ll{CKM)^l'j^ (th)) X 10~^) [29], with the precision of the comparison now being dominated by 
theoretical uncertainties. 

This sample is substantial enough to move to the next step: measurement of polarization am- 
plitudes. This measurement is currently ongoing, and a resolution of 10% is expected. This is 
particularly interesting in the light of current data on other penguin-dominated VV decays, show- 
ing much smaller values for the longitudinal fraction fi than expected in the SM ("polarization 
puzzle") [37]. 

4. Prospects 

At this time, each of the Tevatron experiments has accumulated about 6fb~^ of data. Sev- 
eral current results on flavor physics are based on much smaller samples (1 or 2fb~^), and many 
channels have not been explored at all; this is particularly true for the large hadronic samples at 
CDF. The Tevatron luminosity has been steadily increasing over the past years up to the present 
date. With the current expectation of running with steady luminosity through year 2011, a further 
doubling of statistics is expected, reaching a sample of 10-12fb~^ With this amount of data and as 
the analysis work progresses, many new results are expected. 

The future sensitivity of the analysis has been estimated by CDF, in the conservative as- 
sumption of no analysis improvements, and no use of constrains coming from additional measure- 
ments liko AsL or other external quantities. Fig. 6 show the probability of observing a 5<7 deviation 
from the SM as a function of the true value of j8^, for 10fb~^ of collected data, for a single, or the 
combination of the two experiments (assumed equivalent to CDF). The probability for 5(7 obser- 
vation is very high if is above 0.4 - 0.5 (the probability of 3<7 evidence is of course much higher 
still). It is interesting to note that several BSM physics model do predict values of in this range; 
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Figure 6: Projected probability for observing at the Tevatron a 5(7 deviation from the SM, as a function of 
the true value of ft 



for instance, 4th generation models have formulated predictions in the range 0.5-0.7 [33]. At the 
same time, the related ^ (p(p mode is also going to be studied with larger statistics. Polarization 
will be measured and eventually CP asymmetry. 

The increase of statistics brings not only an improvement in the precision of current measure- 
ments, but also a significant broadening of the range of measurements that are accessible. The 
oscillations were first observed in 2006 by combining several exclusive and inclusive B^ modes, 
and all available flavor tagging algorithms. Today we can observe a clear oscillation signal in a 
single fully reconstructed channel, and using a single tagger. Fig. 7 show the amplitude scan in 
2.8fb-i for 5^ ^ D.TT, using only the same-side kaon tagger (SSKT), exhibiting a clear signal 
at the frequency measured in the original observation [36]. This signal confirms that the tagging 
capability of CDF has not been hampered by the increased luminosity, and it will allow to cali- 
brate tagger dilution from data, rather than having to rely on Monte Carlo. At the same time, a 
more powerful flavor tagging algorithm is being developed, using an ANN to combine information 
from both the same and the opposite side of the B^ in a single tagger, to achieve the best possible 
performance. As a consequence, we are now entering the era of time-dependent CP-asymmetry 
measurements in the B^ system. Examples are the B^ K^K~ mode, and the determination of 
angle y from B^ DsK. 

Other determinations of y will also become important. CDF has performed a pilot measure- 
ment of CP asymmetry in the mode B^ D^K^, where the decays in one of the CP eigenstate 
modes kk and KK ("GLW method"). Using a sample of just lfb~^ it provided a result 
already of comparable resolution with existing results on the same mode from e^e~ B-factories. 
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Figure 7: Amplitude scan for 5^ flavor osciUation in the 5^ DsTt mode, based on Same-Side Kaon 
Tagging. 

Reconstruction of D^cs^ mode was also achieved. From these results, one can extrapo- 

late that the CDF sample with 10fb~^ will be roughly equivalent to 3 ab~^ of e^e~ data, which 
implies it will be the most important measurement at that time, adding on top of the CDF-unique 
measurements in the DsK mode. 

4.1 Charm physics 

A field that is becoming very important at the Tevatron is charm physics. CDF has a huge 
sample of decays from its hadronic trigger, and is accumulating further data at a high rate (as an 
example, the rate of reconstructed Kn events is 4M/year, a rate an order of magnitude larger 
than the next best experiment. Belle). From these data, measurements of both mixing and CPV 
have been performed in the past from small samples. The CDF measurement of oscillations 
in the Kn channel with a 1.5fb~^ sample has produced very similar results to BaBar, both for 
resolution and values measured for the mixing parameters. The current CDF sample should allow 
a 5(7 observation to be performed in a single experiment. 

The measurement of CP violation in kk and KK modes published by CDF in 

2005 using a sample of just 0.13 fb~^ has been the dominant measurement at the time, with a 
resolution of 1.3%, and it is still a significant contribution to today's world average, which has a 
resolution of 0.4% [1]. The analysis in progress on the current data sample is expected to improve 
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the world-average precision by a factor of 2, and reach an ultimate precision at the end of the run 
in the 10~^ range, a very interesting region for possible new physics effects in the charm system. 
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